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Abstract. We present an optimization of the hadron focusing system for a low-energy high intensity con- 
ventional neutrino beam (Super-Beam) proposed on the basis of the HP-SPh at CERN with a beam power 
of 4 MW and an energy of 4.5 GeV. The far detector would be a 440 kton Water Cherenkov detector 
(MEMPHYS) located at a baseline of 130 km in the Frejus site. The neutrino fluxes simulation relies on 
a new GEANT4 based simulation coupled with an optimization algorithm based on the maximization of 
the sensitivity limit on the 6*13 mixing angle. A new configuration adopting a multiple horn system with 
solid targets is proposed which improves the sensitivity to #13 and the CP violating phase 8c p- 

PACS. 14.60.Pq Neutrino mass and mixing, Super Beam, MEMPHYS, SPL. 



1 Physics goals 

The determination of the — > v e oscillation probability 
is one of the main goals of current research in neutrino 
physics. This process contains information on currently 
unknown fundamental parameters of the PMNS matrix 
which describes the mixing between the neutrino mass 
and flavor eigenstates: the #13 angle and the 5c p phase. 
Moreover the ordering scheme of neutrino masses (the so 
called mass hierarchy) could be determined through effects 
induced by neutrino interactions with matter through the 
MSW effect. Given a non null #13, CP violating effects 
generated by Sep 7^ or n would induce a difference in 
the oscillating behavior of neutrinos and anti-neutrinos: 
V(vii —> v e ) 7^ "P(Pii ~ * Ve). The presence of such effects 
involving light neutrinos would support the leptogenesis 
hypothesis as an explanation for the observed matter/anti- 
matter asymmetry in the Universe. 

With respect to — > v T transitions for which the 
mixing is maximal, v e appearance is a sub-leading effect 
suppressed by the smallness of the #13 angle. This imposes 
the use of a very intense and pure neutrino beam and a 
Mton scale low density detector capable of measuring the 
appearance of v e . 

In this work we present a development of the original 
proposal for a high-intensity conventional neutrino beam 
(Super-Beam) produced at CERN and aimed towards a 
440 kton Water Cherenkov detector at the LSM laborato- 
ries in the Frejus tunnel, at a distance of 130 km. 

This relatively short baseline is beneficial in terms of 
suppression of the neutrino flux. Furthermore the sub-GeV 
spectrum is well matched to the region in which quasi- 
elastic interactions are dominant. Water Cherenkov de- 
tectors provide excellent reconstruction for this topology 
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with good rejection capability for the ir° mesons produced 
in neutral current interactions. Furthermore a clean mea- 
surement of CP violation (CPV) is possible since there is 
almost no interplay with effects related to propagation in 
matter. 

First studies [TJ[2|3] were assuming a 2.2 GeV proton 
beam and a liquid mercury jet target associated with a 
single conic horn with a pulsed current of 300 kA. Later 
it was proposed [4] to supplement the system with an 
auxiliary horn (called reflector) enclosing concentrically 
the first horn and operated at 600 kA in order to focus 
also pions produced at larger angles. This scheme was 
adopted in [5] and the horn shape re-optimized using the 
method described in [SJ- Further, the decay tunnel was re- 
optimized using different primary beam energies from 2.2 
up to 8 GeV. Based on the neutrino fluxes of [5] and an 
improved parametrization of the far detector, the physics 
performances of the project were presented in [7[ assuming 
a 3.5 GeV proton kinetic energy. 

With respect to previous studies on this subject we 
propose a new design based on the use of a solid target 
and a single magnetic horn operated with a lower value of 
the pulsed current (300-350 kA). Such a setup simplifies 
the engineering complexity of the system avoiding difficult 
issues such as the containment of the mercury jet in a 
magnetic field free region [8], the challenge of a power 
supply operating at 600 kA and the constraints related to 
mechanical stresses on the horn-reflector system induced 
by the high frequency current pulsing. The capability for 
a solid target to sustain for a reasonable time a power of 
4 MW has not yet been demonstrated. For this reason we 
propose a setup based on a battery of four target-horn 
systems operated sequentially [9]. 

In Se ction [2] we descr ibe the proton sour ce and far 
detector (|2.1[), the target (2.2) and horn (2.3) setup, the 
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method used to calculate the neutrino fluxe s (2.4) and the 
sensitivity to the physical parameters (2.5). Sect ion [3] de- 
scribes a systematic procedure of optimization of the horn 
and decay tunnel based on the sensitivity on sin 2 2r3 13 . 
The neutrino fluxes obtained with an optimized setup are 
the content of Section [4] together with the sensitivities 
to sin 2 2f9 13 and Sep for which we also take into account 
the uncertainty related to hadroproduction correcting the 
simulation with experimental data. Finally Section [5] de- 
scribes two alternative options for the target design. 



2 Setup 

2.1 Proton driver and far detector 

The proton driver we consider is the high power version of 
the Super Conducting Proton Linac (HP-SPL) at CERN. 
The current design study [TU] foresees a proton kinetic en- 
ergy of 4.5 GeV and a beam power of 4 MW at 50 Hz repe- 
tition frequency with a pulse duration of about 400 /is. An 
accumulator ring divides the protons into pulses of ~ 1 /xs 
with a sub-structure of either 6, 3 or 1 bunches per pulse 
pTj . It should be noted that this machine supplemented 
by a compressor ring could constitute the proton driver 
for a future neutrino-factory. 

The far detector which we consider is a water-Cherenkov 
with a 440 kton fiducial mass (MEMPHYS) |Hj. A design 
with a mass of 500 kton is also being considered. Besides 
its role in neutrino oscillation physics, the programme of 
such a detector includes the study of proton decay, atmo- 
spheric and SuperNovae neutrinos. 



2.2 Target 

Due to the low energy of the proton beam pions are emit- 
ted at relatively large angles. This constraint forces to 
place the target inside the horn in order to preserve a 
good collection efficiency. 

In previous studies the liquid mercury target has been 
schematically modeled as a cylinder with 0.75 cm radius 
and 30 cm length corresponding to about two hadronic in- 
teraction lengths (A/). In this study graphite was chosen 
as an alternative material since solid targets composed of 
graphite constitute a proven technology in existing neu- 
trino beams (i.e. T2K and CNGS). It must however be 
noted that other low atomic number (low-Z) materials 
such as Beryllium, Aluminum or AlBeMet® (38% Al, 
62% Be) give similar particle multiplicities and spectra. 
A granular Titanium structure has also been considered 
and is discussed in Sect.O 

In [T3] a comparison of the mercury and the carbon 
target in terms of energy deposition and secondary parti- 
cle yields is drawn. The power deposited in the target at 
4.5 GeV and 4 MW power is about 220 kW and 700 kW 
for graphite and mercury respectively. Pion yields are sim- 
ilar while the most remarkable difference between the two 
targets is the neutron yield which is reduced by a factor 15 
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Fig. 1. Particle yields versus the proton kinetic energy for a 
30 cm long, 1.5 cm diameter mercury target and a 78 cm long 
1.5 cm diameter graphite target with a 4 MW beam power 
according to the FLUKA simulation. 



for the graphite target (see Fig. [I]). This reduction is ben- 
eficial in view of the radiation damage caused by neutrons 
on the horn|14j. In this work we consider a cylindrical tar- 
get with a radius of 1.5 cm, a length of 78 cm and a density 
of 1.85 g/cm 3 corresponding to roughly 2 A/. 



2.3 Horn 

In the previous study the horn shape is designed in order 
to focus pions of ~ 600 MeV to produce neutrinos with 
energies at the oscillation maximum (E v ~ 260 MeV). As- 
suming a point-like source and a cylindrical surface close 
to the target the optimal shape is roughly conical. 

In [13] we have shown that the conical horn optimized 
for the liquid mercury target is not suited to a long solid 
target since the sin 2 2f5i3 sensitivity limit for Sep < tt 
degrades significantly. The effect was found to be caused 
by the contamination of v e in the beam from cascade 
decays of positive pions {tt + — > fi + — > e + v e D^) and kaons 



(K+ 



which are not effectively defocused when 



they are produced in the forward direction^] 

A new horn model inspired by the one used for the 
MiniBooNE beam [IS], having a larger acceptance for for- 
ward produced pions, proved to be quite effective with 
respect to the need for a reduced contamination from 
wrong-charge pions. The generic layout of the horn which 
we will use in the following (forward-closed horn here- 
after) is shown in Fig. [2] The geometry has been pa- 
rametrized in terms of the longitudinal displacement of 
the target with respect to the horn (zq S ), the lengths of 
the longitudinal (£1,2,3,4,5) an d radial (-Ri,2,3) sections, 
two curvature radii (r^a) and the conductor thicknesses 
{h, 2,3,4)- 

We consider two scenarios for the horn/target system: 



1 The equivalent effect is weaker in the v beam due to the 
combined effect of the reduced cross section of v e and the fact 
that negatively charged hadrons are less abundantly produced. 
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probability for the decay v to be emitted through a 
surface of area A at a distance L is: 



Fig. 2. Parametrization of the forward-closed horn geome- 
try. The system has a rotational symmetry around the z axis. 
Figure shows a cut in the r — z plane. 



— Integrated conducting target. The target could be an 
integral part of the inner horn conductor allowing for 
higher magnetic fields at the target exit (7?i +t\ = 
R*b = 1.5 cm). Such a setup was studied in [16] taking 
into account energy deposition in the target due to the 
primary beam and the resistive (Joule) losses caused 
by the horn current which passes the common target- 
inner conductor region. 

— A target separated from the horn. We will present how 
the physics performances vary when the horn inner ra- 
dius is increased. In this case the cooling system could 
be realized with a helium based annular duct flow cool- 
ing as done in T2K [17] . 

Taking advantage of the small transversal dimensions 
of the horn, the idea of using a battery of four horns in 
parallel has been proposed. This arrangement would imply 
reduced stress on the targets via lower frequency pulsing 
(12.5 Hz). This choice would bring the beam power on 
each target to 1 MW which is currently considered as a 
viable upper limit for solid targets operations. 



2.4 Neutrino flux simulation 

The neutrino energy spectra are calculated using a proba- 
bilistic approach in order to obtain reliable results in a rea- 
sonable amount of time using samples of ~ 10 6 simulated 
protons. The probability that the neutrino will reach the 
far detector is calculated at each particle decay yielding 
neutrinos with analytic formulas [BIITI ITBIH^] . The proba- 
bility is then used as a weight factor in the calculation of 
the neutrino energy spectrum. The neutrino sources con- 
sidered in the simulation are: 

1. Two body decays of charged pions and kaons: 

{n + ,K + } —> and charged conjugate (c.c.) with£ = 
/i, e. Given a parent traveling in the laboratory frame 
with an angle a with respect to the beam axis, the 



V 2 



A 



AttL 2 [p cos a - If 



where f3 and 7 represent the relativistic factors of the 
parent it or K in the laboratory frame (LF). In this 
approximation, the decay region is seen as point-like 
from the far detector at the distance L. 
2. Muon decays: fi + — > e + v e v^ + c.c. The differential 
probability for the and v e to be emitted parallel to 
the beam axis with an energy E v in the LF is: 



dv = 2(j -n^f 1 cose*) 

dE v m P 7 M (l + /3 M cos0*) 



where 9* is the angle between the v and the [i direc- 
tions in the [i rest frame (RF). Denoting the v energy 
in the ,u RF as E*, x = 27?* /m M and f = 2x 2 (3 - 2x), 
h = 2x 2 (l - 2x) for the / = /1 = I2x 2 (l - x) 
for the v e . 77^ is the muon longitudinal polarization: 



77,^ = ,/l - ( :> ^sin0*) 2 where 0* is the angle of the 

\x with respect to the beam axis in the parent (p = tt 
or K) RF. The sign in front of 77^ becomes + for 
decays. This differential probability is integrated 
in E v sampling the energy which is available for the v 
uniformly in bins of 1 MeV width. In order to reduce 
the statistical fluctuation which would arise from the 
limited number of muons decaying in the decay tun- 
nel (/37CT ~ 6.3 km for = 1 GeV/c) each time a /i 
is produced the probability of its decay is introduced 
as a weight. This is calculated from the length trav- 
eled inside the decay tunnel assuming a straight line 
propagation. 

3. Kaons semi-leptonic three body decays: K — > nlvi 
with I = {/i,e} and K = {K ± ,K }. The probabil- 
ity to reach the detector is computed as: 



dV 
dE„ 



(m K — m w - to;)7a'(1 + Pk cos0*) 



V 2 (K) 



6* is the angle between the v and the K directions in 
the K RF. This differential probability is integrated 
in E v sampling the energy which is available for the v 
uniformly with a fixed number of 200 points. f{E*) is 
a parametrization of the distribution of the energy of 
the v in the K RF which was derived from [2"Q] . 

The distribution of the secondaries at target exit ob- 
tained with the FLUKA [5T] generator is used as an ex- 
ternal input to a GEANT4 [22 simulation derived from a 
GEANT3 code developed in [B]. The target, the horn with 
its magnetic field and the decay tunnel are fully simulated 
within GEANT4. Alternatively GEANT4 can be used to 
simulate also the interactions of primary protons in the 
target: this option was used as cross check (Sect. |4~2| ). In 
order to cross-check and validate the new GEANT4-based 
software, a comparison has been done with the fluxes ob- 
tained with GEANT3. The fluxes obtained in the two 
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frameworks are in good agreement both in terms of nor- 
malization and shape |13j . Further cross-checks included 
the correct implementation of the decay branching ratios, 
a comparison with an independent code and a check based 
on direct scoring of the emitted neutrinos. 

2.5 Sensitivities 

The sensitivities for the measurement of the oscillation 
parameters $13 and Sep are obtained with the help of 
GLoBES 3.0.14 [53]. The parametrization of the MEM- 
PHYS Water Cherenkov detector [T5] at the level of physics 
performance (efficiencies, background rejection, etc.) is 
implemented in the public AEDL file SPL.glb which is 
distributed with GLoBES. The results derive from the 
analysis algorithms developed in the context of the Su- 
perKamiokande detector. The efficiency for the reconstruc- 
tion of v e events is 70%, the fraction of misidentified muon 
events is 0.054% of the v,, events and the neutral cur- 
rent ir° background is 0.065% of the un-oscillated v^ c 
and 0.25% of the un-oscillated v^ c . A detailed descrip- 
tion can be found in [TJ. In the following we assume a 
sharing between neutrino and anti-neutrino running of 2 
and 8 years, respectively, as it was done in [TJ. A system- 
atic error of 5% on the fluxes is assumed in the calculation 
of the limits. 

The sensitivity limit is defined by generating the event 
rates in the assumption of a null #13 and performing a fit 
with finite values of #13. In particular, in the optimization 
process (Sect. [3]), we assumed the normal hierarchy and 
the A\ 2 was calculated neglecting the correlations of the 
parameters with the glbCliiSys routine (see |23jV 

3 Sensitivity-based optimization 

The approach which was followed in the optimization of 
the forward-closed horn and the decay tunnel uses the fi- 
nal sin 2 2#i3 sensitivity, i.e. the final physics performance 
as a guiding principle in the ranking of the configurations 
under scrutinjj^] In the evaluation of this quantity a com- 
plex set of relevant factors are given as an input: the nor- 
malization and shape of each neutrino flavor, the running 
time in the positive and negative focusing mode, the en- 
ergy dependence of the cross sections, the backgrounds in 
the far detector and its response in terms of efficiency and 
resolution. 

Given the well known dependence on the sin 2 29^ limit 
on the Sep phase, we introduced the quantity A defined as 
the <5cp-averaged 99% C.L. sensitivity limit on sin 2 2#i 3 
in units of 10~ 3 : 

10 3 f 2n 

A= — / \ 99 {S C p)dS C p (1) 
27r Jo 

This quantity has been used as a practical way of defining 
with a single number the quality of the focusing system. 

2 The possibility to use the sensitivity to CP violation as a 
guiding line will be the object of a future study. 



The key parameters defining the horn and tunnel ge- 
ometry are randomly sampled within specified ranges and 
the correlations with the figure of merit A studied. 



3.1 Beam composition and achievable limits 

We started by studying the correlation between the rel- 
ative normalizations of and v e fluxes and the corre- 
sponding variations of A. For this exercise we chose a spe- 
cific horn configuration providing a A ~ 1 and scaled the 
v ai Vu fluxes globally by a factor c M with 0.5 < c M < 2. 
The same scaling was applied to the v e and v e fluxes with 
an independent factor c e within the same limits. 

The parameter A improves (i.e. decreases) by a factor 
of about 3 when the v e is reduced by a factor 2 and 
increased by the same factor while it worsens by a factor 
2.4 when the variation of the fluxes is done in the other 
sense. The iso-sensitivity levels follow quite closely a law 
c n = -Jce- This behavior reflects the fact that the exper- 
iment significance S / vB in terms of observed events in 
the far detector, to which A is related, is approximately 
invariant for a scaling of the fluxes (c„, c e ) in which c„ = 

We have studied the correlation between the integral 
fluxes and v e with the full simulation under variations 
of the decay tunnel geometry only with the horn configura- 
tion which will be described in Sect. 13.21 The correlation is 
strong since the bulk of the v e are generated by the same 
decay chain producing (n + — > ^ + ^, /i + — > e + v e v^ 
+ c.c). The dependence follows roughly a quadratic law: 
<j>{v e ) ~ 0(f M ) Q , a ~ 2. Deviations occur especially at 
small fluxes (<M*v) - 3.5- 10 14 i//100m 2 /year at 100 km) 
below which a < 2. This sets approximately the threshold 
above which increasing the flux by tuning the decay 
tunnel does not improve the sensitivity. 



3.2 The optimization procedure 

The parameters of the forward-closed horn and of the 
decay tunnel were sampled with uniform probability dis- 
tributions imposing the configuration to be geometrically 
consistent ( "iteration- 1"). The decay tunnel length (L tun ) 
and radius (i? tun ) were restricted in the intervals [35, 45] m 
and [1.8,2.2] m respectively. The maximal length and ra- 
dius of the horn were limited to 2.5 m and 80 cm in order to 
maintain a compact design which allows to use a battery of 
four target-horn stations in parallel. Moreover, the inner 
radius R\ was limited in [1.2, 4] cm, the lower limit cor- 
responding to the "integrated target" limit. Further con- 
straints on the parameters were L\ > 50 cm, L5 < 15 cm 
and —30 < Zg S < cm. 

The stability of the value of A related to statistical fluc- 
tuations (2 • 10 5 tracks per configuration) was estimated 



3 The signal (S), the v e appearance, is a fixed fraction of 
the flux and the background (B) is mostly generated by the 
intrinsic v R component in the beam. 
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Fig. 3. Distribution of A for iteration-1 (continuous, 3000 
configurations), iteration-2 (dashed 1000 configurations), and 
iteration-3b (dotted, 1000 configurations). See the text for the 
definition of the samples. 
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Fig. 4. Distribution of the horn geometrical parameters Ri 
(inner radius) and Ri + R2 (alcove radius) for the inclusive 
sample (empty histogram) and for the S subsample (filled his- 
togram). See text. 



by repeating the simulation using an identical set of pa- 
rameters with independent track samples. We find a r.m.s. 
of the integral and v e fluxes respectively at the level of 
3% and 5% of the central value. 

The distribution of A is shown in the continuous his- 
togram of Fig. [3] Studying the distributions of the in- 
tegral fluxes for positive and negative focusing for the 
sub-sample (S) of configurations providing good limits 
(A < 1.05) we observe an enrichment in configurations 
yielding high fluxes with a low v e and contamina- 
tion as expected. Small v e contamination in the 9u flux 



225 cm! 



56 cm' 



is confirmed to be also important (see Sect. 2.3). On the 
other hand the mean and the root mean square of the 
energy spectrum do not show a visible correlation with A. 

The statistical distributions of the geometrical param- 
eters for the inclusive sample and the S subsample were 
compared to obtain an indication of the parameters which 
are more effective in increasing the sensitivity. Despite the 
averaging on a large set of parameters variations, the dis- 
tributions show interesting features for the following pa- 
rameters: the horn inner radius (i?i), which has a strong 
preference for small values; the target position (zq S ), for 
which very low values are disfavored; the radius of the 
"alcove" (Ri + R2) whose optimal values cluster at about 
20 cm and the decay tunnel length (L tun ) which tends to 
give better performance for small values. Two of these pa- 
rameters and R± + R2) are shown as an example in 
Fig. [4} 

According to these indications a second scan ( "iteration- 
2") was performed after fixing the parameter i?i at 1.2 cm 
and restricting the ranges of variation: 20 < i?i + R2 < 
22 cm, 30 < L tun < 40 m, -15 < z* s < cm. The dis- 
tribution of A for this new sample is shown by the dashed 
histogram of Fig. [3j With respect to the same distribution 
for the previous scan (continuous line) a shift towards bet- 
ter limits by about 20% is achieved. The horn shape for 
the configuration giving the minimum value for A in the 
iteration-2 is shown in Fig. [5| 




Fig. 5. Shape of the optimal horn. 
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Table 1. Parameters of the optimized system expressed in cm. 



Finally we kept the horn shape fixed (Fig. [5| and per- 
formed a further tuning of the decay tunnel length and 
radius. The scan was done in two regions to better con- 
strain the position of the minimum: 0.5 < R tun < 2.5 m 
and 10 < L tun < 60 m ( "iteration-3a" ) and 1.5 < R tun < 
4.5 m and 15 < L tua < 35 m ("iteration- 3b"). The dis- 
tributions of the parameter A for the latter sample, which 
lies in the neighborhood of the minimum, is shown by the 
dotted histogram of Fig. [3] With respect to the initial dis- 
tribution an improvement of 25-30% is obtained. 

The dependence of A on the decay tunnel variables 
can be reasonably fitted with a quadratic function: A = 



0.94+2. M0- 4 (L tun [r 



-31. 



-2.4-10- 2 (i? tun [m]-2.9)' 



(Fig. [6j left). Since the minimum is relatively broad we 
chose L tun = 25 m and R tun =2mas central values based 
on practical considerations related to the excavation and 
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shielding of large volumes. This compares to the previous 
values of 40 m of length and 2 m of radius. 




tunnel length (m) current (kA) 



Fig. 6. Dependence of A in the L tml - R tun plane (left) and 
in the (Ri, i) plane (right). 

We have also observed that an increase in the current 
(between 300 and 400 kA) tends to systematically produce 
better sensitivity limits. The interplay between the current 
and the inner horn radius in terms of the quantity A, after 
having fixed the optimal horn and decay tunnel has been 
studied. Data are well fitted by a linear function in two di- 
mensions: A = (9.2- 0.81 -/[100 kA])/(7.3-0.37-Ri[cm]) 
(right-hand plot of Fig. [6]). The effect of increasing the 
current, i.e. a stronger magnetic field in the vicinity of 
the target, is physically equivalent to decreasing the mini- 
mum horn inner radius. In this way, at first approximation, 
working at constant I/Ri (oc B), allows to stay at fixed 
sensitivity. 

In the following we will take the configuration with a 
current of 300 kA and integrated target as the baseline 
choice. It can be noticed however that, in the case of a 
non-integrated target, if the minimum horn radius would 
need to be extended up to 4 cm to accomodate for the 
cooling system, similar performances could be recovered 
by increasing the current up to about 400 kA (right-hand 
plot of Fig. [6|. 

4 Results with the optimized setup 

The distribution of the longitudinal coordinate (z) of the 
decays in flight of charged pions and kaons with the opti- 
mized configuration in the positive focusing mode is shown 
in Fig. [7J The continuous histogram on the left-hand side 
represents the z distribution of the exit point of positive 
pions from the target. The focusing effect is visible in the 
different shape observed for positive and negative mesons. 
The relative yield and lifetime of pions and kaons can be 
also appreciated. 

The focusing effectiveness of the system can be char- 
acterized by comparing the and spectra at the far 
detector obtained with null current and the nominal one. 
Working in the positive focusing mode the flux is in- 
creased of a factor 6.5 and the flux decreased by a fac- 
tor 5.4 with respect to the situation in which no focusing 
would be applied (Fig. [8]) . 

The obtained v e and charged conjugate (c.c.) neu- 
trino fluxes are shown in Fig. [9] for positive (left) and neg- 
ative focusing (right) runs. They correspond to 5.6 • 10 22 
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Fig. 7. Distributions of the longitudinal coordinate (z) of the 
decays in flight of positive and negative pions and kaons with 
the optimized configuration. 
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Fig. 8. Effect of the focusing, and P M spectra in the positive 
focusing mode with horn current off (left) and on (right). 



protons on target (p.o.t.)/year (4 MW • 10 7 s at 4.5 GeV) 
and are calculated at a reference distance of 100 km over a 
surface of 100 m 2 . The fractions of v^, v^, v e and v e with 
respect to the total are (98.0%, 1.6%, 0.42%, 0.015%) and 
(95.3%, 4.4%, 0.28%, 0.05%) for the positive and negative 
focusing modes respectively. In positive (negative) focus- 
ing mode the v e (p e ) fluxes are dominated by muon de- 
cays: 82% (90%). The c.c. fluxes receive instead a large 
contribution from kaon three body decays (81% and 75% 
in positive and negative focusing respectively) with muon 
decays from the decay chain of "wronge charge" pions at 
low energy contributing for the rest. The fluxes are pub- 
licly available [24] . 

In Fig. [lO] the fluxes obtained with the optimized horn 
are compared to those obtained with the original double 
conical horn with currents of 300 and 600 kA associated 
with a mercury target and published in [7J. The za, and 
v e energy spectra are shifted to higher energies with an 
increase in statistics particularly around 500 MeV. The 
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Fig. 9. Neutrino fluxes obtained with the optimized horn and 
decay tunnel in positive (left) and negative (right) focusing 
mode. The table contains the integral neutrino flux per year for 
each flavor at a distance of 100 km over a surface of 100 m 2 . The 
fluxes were obtained with a sample of 10 7 simulated proton- 
target interactions. 
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Fig. 10. Comparison of the neutrino fluxes obtained with the 
new design (continuous line) and the previous one [7] (dotted 
line) . 
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Fig. 11. Variation of the ratio of fluxes obtained with a set 
of four displaced horns (0(r)) with respect to a single centered 
horn (<j>(0)) as a function of the displacement r. Starting from 
the left the first line represents the horn radius, the second and 
third lines the minimum and maximum possible displacement. 
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Fig. 12. Polar angle versus momentum distributions for pos- 
itive pions at target exit (left) and after reweighting with the 
muon neutrino spectrum (right). The grids represent the bin- 
ning used in the HARP analysis. 



only 10%. The baseline configuration with horns as central 
as possible (i.e. in mutual contact, r ~ R horn V2) causes a 
minor loss of of the order of 1-2%. The presence of a 
magnetic field in all the horns simultaneously or in each 
horn separately does not alter significantly the fluxes. 



flux is enhanced also in the proximity of the oscillation 
maximum at 260 MeV where the v e flux is reduced by a 
similar fraction. The wrong-CP component {p e , v^) on the 
other hand is reduced by more than a factor two. 



4.1 Parallel horns 



Fig. 11 shows dependence of the flux obtained using 
four horns as a function of their radial displacement (r) 
with respect to the decay tunnel axis. Even in the most 
pessimistic case where the horns should be placed at the 
decay tunnel edge (r = R tun — i? horn ), the flux reduction is 



4.2 Hadroproduction uncertainties 

In order to characterize the subsample of produced pions 
giving the most important contribution to the neutrino 
flux, the distribution of the momentum (p) and polar an- 
gle (9) of positive pions was considered in bins of the neu- 
trino energy. The pion distributions in each neutrino en- 
ergy bin were averaged taking the flux in the correspond- 
ing bin as a weight. In Fig. 12 the p versus 9 distribution 
of positive pions at target exit is shown before (left) and 
after (right) the above-mentioned reweighting. While the 
inclusive sample peaks at p ~ 0.5 GeV/c and 9 ~ 0.5 
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rad, the pions which mainly contribute to the flux peak 
at higher momentum (p ~ 1.2 GeV/c) and smaller angle 
(9 ~ 0.2 rad). The grids superimposed to the distributions 
represent the binning used in the analysis of the HARP 
hadro-production experiment. The interesting phase space 
region is mainly covered by the so called "small angle" 
analysis for which only "thin" target]^] data exist at 5 GeV 
proton kinetic energy for carbon 25J. 

We used these data to assess the size of the systematic 
effects introduced by the uncertainty on hadro-production. 
HARP pion differential cross sections have been compared 
with the expectations from GEANT4-QGSP and FLUKA. 
In the comparison the geometry of the targets used in the 
HARP experiment has been modeled with the two simu- 
lation program^] 



The comparison is shown in Fig. 13 for GEANT4- 
QGSP only. The absolute normalization and the trend 
of the experimental data is not too far from the simula- 
tion even though large discrepancies are observed in sev- 
eral bins. Overall the QGSP model tends to be slightly 
closer to the data predicting softer momentum spectra. 
As a systematic check we re-calculated the fluxes: 1) us- 
ing the GEANT4-QGSP model in place of FLUKA; 2) 
using FLUKA with reweighting factors derived from the 
comparison of FLUKA and HARP data. The reweighting 
factors in the considered momentum and polar angle bins 
range between 0.14 and 1.6. The reweighted spectrum 
is peaked at lower energies where it exceeds the original 
one while the high energy tail is reduced. We will show the 
difference of these assumptions on the primary pion spec- 
tra at the level of sensitivities in the following section. It 
must be noted that with this procedure we are neglecting 
possible effects related to the use of a long target and the 
effect of the slightly different proton kinetic energy. 



4.3 Discovery potential for 6>i 3 ^ and CP violation. 

Event rates in the MEMPHYS detector are shown in Fig. [14] 
for sin 2 2^i3 = 0.01 and 8cp = 0. In positive focusing 
mode the main background comes from the intrinsic beam 
contamination (84%) followed by muon events misidenti- 
fied as electron events (9%) and neutral current interac- 
tions (NC, 7%). In the negative focusing mode the intrinsic 
beam contamination amounts to 65% (35% v^ c and 30% 
Ve C ), the misidcntificd event to 6% and the NC interac- 
tions to 28%. 

In Fig. [15] the curves define the region in the space of 
true (sin 2 2#i 3 , Sep) for which a discovery of sin 2 29 13 7^ 
is possible at more than 3 a. The limit as a function of 
8c p is defined as the true value of sin 2 2#i3 for which a fit 
assuming #13 = gives a Z\% 2 = 9 as done, for example, 



4 A cylinder with 1.95 cm length and 1.5 cm radius. 

5 The formula which we used reflects the approach used in 



d (T 



■jfeNij where JVy 



arc 



the HARP publications: d > ,,, v 

the number of events in the bin (i, j) at true level, N pot is the 
number of simulated protons interactions in the sample, p the 
density of the graphite target, A the atomic number of carbon 
and t the target length. 
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Fig. 14. Event rates in the MEMPHYS detector for 
sin 2 26*13 = 0.01 and Sep = for 8+2 years of v+v running. 



in [7] . In the fit a minimization on all the other oscillation 
parameters is done and the two hypotheses for the mass 
hierarchy are taken into account. 

The limit obtained with the previous setup associated 
with the mercury target is shown by the dash-dotted curve 
while the new limits are represented as a hatched band. 
The upper edge of the band (continuous line) refers to the 
FLUKA model of hadro-production, the lower edge (dot- 
ted) to the GEANT4-QGSP model, the one lying (mostly) 
in the middle (long dash-dotted) is obtained after reweight- 
ing FLUKA with the HARP data. 

The CPV discovery potential at the 3 a level is shown 
in Fig. |16| discovery is possible in the region defined by 
the couple of curves. This means that in that region of 
the true (sin 2 20i 3 , Sep) plane a fit done under the CP 
conserving hypotheses (5c p = 0,7r) gives for both choices 
a A\ 2 > 9. The meaning of the symbols is the same as 
for the 9 13 plot. The new limits generally improve those 
obtained with the previous design both for 9i 3 and CPV 
discovery. Considering the results obtained with FLUKA 
without reweighting (solid line) in the region dominated 
by the v run (6 < n) the best limit is improved by about 
4% (from 1.35 to 1.30 10 -3 ) and in the complementary 
^-driven region by about 25% (from 1.08 to 0.80 10~ 3 ). 
The same criterion for the CPV discovery gives and im- 
provement of 10% (from 1.50 to 1.35 10~ 3 ) and 25% (from 
1.65 to 1.23 10~ 3 ). 



5 Alternative target designs 

As undergoing studies in the context of the EUROnu de- 
sign study have shown technical challenges for the inte- 
grated horn-target system (IT), we studied the perfor- 
mances of two additional configurations assuming: 

— a graphite target separated from the horn (ST); 

— a granular target composed of titanium spheres with 
diameters of O(mm) (PB, packed pebble bed target). 

while keeping the target geometry unchanged. 
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Fig. 13. The HARP data points [55] (bullets) are compared to the GEANT4-QGSP model (solid red histograms) for ir + (upper 
row) and n~ (lower row). Data correspond to 5 GeV protons impinging on a thin carbon target. 
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Fig. 15. #i3 discovery potential at 3 a level. We set the "true" 
values of the parameters at Am\ 2 = 7 ■ 10~ 5 eV 2 Am^ — 
+2.43 ■ 10" 3 eV 2 (normal hierarchy), 6 12 = 0.591 and 6> 23 = 
7r/4. We included these parameters in the fit assuming a prior 
knowledge with an accuracy of 10% for #12, 623, 5% for Zim§i 
and 3% for at Am\ 2 at 1 a level. The running time is (2v+8i>) 
years. 



In the two cases we set for the inner radius of the horn 
(i?i) a value of 3 cm. The gap of 1.5 cm between the tar- 
get and the horn is intended to accomodate the infrastruc- 
ture needed for the target cooling. As we showed earlier 
(right hand plot of Fig. [6| the loss in performance due 
to an increase in the horn inner conductor radius can be 
compensated by an increase in the pulsed current. In the 
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Fig. 16. CP violation discovery potential at 3 a level. See also 
Fig. [15] 

following we will assume a value of 350 kA. This current 
is not far from the one used in T2K horns which amounts 
to 320 kA. 

The granular target design was originally proposed in 
26 . Thanks to the favorable surface to volume ratio and 
the possibility to flow transversely the coolant within the 
interstices of the spheres (i.e. a high pressure flow of He 
gas), this target is expected to have a good behaviour even 
under the extreme irradiation conditions we are dealing 
with (1 MW beam power). 

The pion yields for the PB target have been calculated 
assuming a density reduction of 74% which is obtainable 




using the optimal packing of the spheres (central pack- 
ing). We have verified that the reduction in density related 
to the granular structure tends to produce a longitudinal 
profile for the emitted pions which resembles the one ob- 
served for the carbon target which has been used in the 
optimization. This effect arises from the density reduction 
and the shorter hadronic interaction length of Titanium 
with respect to graphite. Even if the achievable density 
reduction fraction would be lower than the assumed value 
the inpact on the sensitivity would not be strong thanks 
to the fact that we are using a target which is relatively 
long in units of A/. 

The discovery potentials for 6*13 ^ and CP violation 
for these two configurations (SB and PB) are compared 
to the ones obtained with the former design based on a 
mercury target (HG) and to the performance with the 
integrated target (IT) in Figs. 17 and 18 With respect 
to the IT design the solution witna separated monolithic 
graphite target and increased current (ST) introduces a 
moderate worsening of the limits; the PB solution gives 
practically unchanged performance for S > tt and some 
improvement for 5 < tt. This effect can be understood in 
terms of a higher symmetry in the production of positive 
and negative pions for Titanium. 

We can conclude that a PB target in association with 
the optimized horn represent possibly the most appealing 
solution in terms of both physics performance and engi- 
neering. 



horn system is foreseen to be replicated four times and 
operated in parallel. With respect to the previous setup we 
abandoned the liquid mercury target option, the 600 kA 
pulsed reflector and shortened the decay tunnel (formerly 
40 m). The horn shape has been adapted to focus pions 
emerging from a long target through a systematic opti- 
mization procedure based on a forward-closed model and 
a ranking of the configurations relying on the sensitivity 
limit on $13. The scenario involving four horns in paral- 
lel has been validated in terms of the achievable neutrino 
flux. The final results of the sensitivities to sin 2 2813 and 
Sep improve those obtained in [6] [7]. The uncertainty on 
hadro-production has been addressed, for a graphite tar- 
get, at the level of sensitivities by exploiting the data of 
the HARP experiment and different models (FLUKA and 
GEANT4-QGSP). 

A strategy considering higher energy proton drivers 
and longer baselines is being addressed with a similar sim- 
ulation and optimization approach and will be the subject 
of a future publication. 

I wish to thank M. Zito with whom I had a fruit- 
ful and intense collaboration, A. Cazes, J. E. Campagne, 
authors of the original GEANT3 simulation, and C. Bo- 
beth, M. Dracos, P. Cupial, C. Densham, M. Mezzetto, 
M. Bonesini. I acknowledge the financial support of the 
European Community under the European Commission 
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6 Summary and conclusions 

In this work the neutrino beam production for the SPL- 
Super Beam from CERN to Frejus has been optimized 
for the measurement of the oscillation parameters #13 and 
5qp aiming, at the same time, to a simplified setup in 
terms of technical feasibility. 

The envisaged scenario adopts as a baseline the use 
of a solid target, a single horn pulsed with a current of 
300 kA and a decay tunnel length of 25 m. The target- 
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